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Abstract 

We study a simple class of unitary renormalization group (RG) transformations governed by a 
parameter / in the range [0,1]. For / = 0, the transformation is one introduced by Wegner in 
condensed matter physics, and for / = 1 it is a simpler transformation that is being used in nuclear 
theory. The transformation with / = diagonalizes the Hamiltonian but in the transformations 
with / near 1 divergent couplings arise as bound state thresholds emerge. To illustrate and diagnose 
this behavior, we numerically study Hamiltonian flows in two simple models with bound states: 
one with asymptotic freedom and a related one with a limit cycle. The / = transformation 
places bound-state eigenvalues on the diagonal at their natural scale, after which the bound states 
decouple from the dynamics at much smaller momentum scales. At the other extreme, the / = 1 
transformation tries to move bound-state eigenvalues to the part of the diagonal corresponding 
to the lowest momentum scales available and inevitably diverges when this scale is taken to zero. 
Intermediate values of / cause intermediate shifts of bound state eigenvalues down the diagonal 
and produce increasingly large coupling constants to do this. In discrete models, there is a critical 
value, f c , below which bound state eigenvalues appear at their natural scale and the entire flow to 
the diagonal is well-behaved. We analyze the shift mechanism analytically in a 3x3 matrix model, 
which displays the essense of this RG behavior and we compute f c for this model. 
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I. INTRODUCTION 



Wilsonian renormalization group transformations typically eliminate (integrate out) de- 
grees of freedom whose energy is much higher than those of interest, replacing them with 
effective scale-dependent interactions. Such transformations allow one to tune an effective 
theory resolution, focusing on essential degrees of freedom and interactions at any scale of 
interest. Anticipated by Kadanoff's block spin transformation [JJ, early transformations ex- 
plicitly reduced the number of degrees of freedom by lowering cutoffs on energy [2j [3] . Glazek 
and Wilson introduced a similarity renormalization group (SRG) procedure which instead 
uses transformations that do not remove any degrees of freedom but eliminate couplings 
between disparate energy scales [U El E] . After regularization and identification of necessary 
counterterms, the SRG procedure eventually produces a renormalized, band-diagonal matrix 
representation of the Hamiltonian. Independently, Wegner introduced non-perturbative dif- 
ferential flow equations that unitarily transform Hamiltonian matrices to the band-diagonal 
form IB], and Wegner's flow equations can be employed in the SRG procedure. 

A simplified version of Wegner's transformation [5] has been successfully applied to a 
number of nuclear few-body problems [TQl [Hj . Nuclear many-body calculations are plagued 
by strong nucleon-nucleon correlations due to a hard core and strong short-range tensor force, 
so perturbative and variational methods converge poorly. The new simple transformation 
produces universal nucleon-nucleon interactions with drastically improved perturbative and 
variational behavior. It has also been shown that one can apply SRG with Wegner's flow (or 
a suitably altered flow) for studying the connection between asymptotic freedom and limit 
cycles [32] and it has been suggested that an infrared limit cycle may exist in QCD |13j . 
These examples indicate that convergence properties of SRG transformations are relevant 
to the theory of particles over a broad range of energies. 

We find that convergence properties of the simple transformation are worse than those 
of Wegner's transformation. The simple transformation tends to diverge whenever the SRG 
parameter approaches the momentum scale at which a bound state is formed in the the- 
ory. This effect may present no problem for applications in low-energy nuclear physics as 
long as the SRG parameter stays larger than the momentum scales at which formation of 
bound states occurs. However, such fortunate conditions may not be available if the Efimov 
effect [TU US] in the three-nucleon problem [T71 HBJ EB [20] cannot be avoided. The 
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Efimov effect shows up as a limit cycle. We will see that Wegner's transformation is capa- 
ble of resolving limit cycle behavior, while the simple transformation diverges as the first 
high-energy cycle is resolved. 

Transform a Hamiltonian H(0) using a unitary operator U(s), 

H(s) = U(s)H(0)U^(s) } (1) 

where s is the SRG flow parameter. We want to choose U(s) so that H is diagonalized as 

s — > oo (band-diagonalized for finite s). We choose s = for the initial value, so H(0) 

can be thought of as an input bare regularized Hamiltonian with all required and properly 

adjusted counterterms (such Hamiltonians are established using the same SRG procedure 

but this aspect is not in our focus here). 

Taking derivatives of both sides of Eq. (1), we see that H(s) evolves according to 

dH{s) 
ds 

with 



[v(s),H(s)}, (2) 



Vis) = ^r^OO = -vHs) ■ (3) 
Choosing r](s) specifies the transformation. We study transformations that mix two simple 
choices, 

V (s) = [D(s),H(s)}, (4) 
where D(s) is the diagonal part of H(s) in momentum representation, and 

v( s ) — [T, H(s)] , (5) 

where T is a fixed matrix, here chosen to be the kinetic energy. If one considers T to be an 
arbitrary Hq that has a known spectrum, the interaction is the remaining part, Hi = H — Hq. 

Using D in rj was Wegner's initial choice and this transformation has been studied ex- 
tensively [8]. The use of T was explored perturbatively [9] and then shown to effectively 
decouple low- and high-momentum scales in a universal characterization of the nucleon- 
nucleon interaction [TU1 [EE] . Universality in effective nucleon-nucleon interactions was dis- 
covered earlier [211 E2] using the same transformation Wilson used in his initial numerical 
RG calculations [2]. 

We will see that using T instead of D produces singularities starting at bound state 
thresholds and limits how far the transformation can be run. But we stress that if it is not 
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run too far, for any phenomenologically tuned nucleon-nucleon interaction H(0) the low- 
energy part of H(s) is nearly universal. This means that all infrared-constrained nucleon- 
nucleon interactions with bare cutoffs at or even well above 500 MeV collapse onto a nearly 
universal infrared Hamiltonian after the SRG parameter, A = 1/ ^/s, playing the role of 
effective cutoff, is evolved to well below the bare cutoff. These evolved potentials disagree 
only on the high-energy part of H(s) , which is not constrained by low-energy physics [TT] . 
Introducing a parameter / which takes values in the range [0, 1], we define 

G f (s) = fT + (l-f)D f (s), (6) 

and write 

±H f ( 8 ) = [F{H f (s)},H f (s)}, (7) 

H f (0) =H, (8) 

where H is an initial Hamiltonian matrix and the SRG generator matrix takes the form 

F{H f (s)} = [G f (s),H f (s]}, (9) 

and thus defines 77 that depends on / and interpolates between the two cases from Eqs. Q 
and (|5| when / varies from to 1. For explicit calculations one uses a basis in which T is 
diagonal since it is in this representation that H(s) is driven towards band-diagonal form. 

The generator F{Hf(s)} is guaranteed to bring the Hamiltonian to diagonal if some 
sufficient and easily verifiable conditions are satisfied. Such conditions will be discussed 
after we introduce all details that count in the derivation, in the context of two examples of 
basic interest in physics. 

Namely, we use the family of generators F{Hf(s)} to evolve Hamiltonian matrices that 
exhibit asymptotic freedom and limit cycle behavior. The asymptotically free matrix model 
is easily derived from the nonrelativstic Schrodinger equation in two dimensions with a delta- 
function potential. Isolate the angular momentum zero states and discretize the momentum 
so that p — > b n p with b > 1; include the appropriate weights to reproduce the momentum 
representation bound-state integral equation in the limits where n is allowed to be any 
integer and b 1. Introduce cutoffs so that M < n < N and you have the matrix we use 
to illustrate asymptotic freedom (M is a large negative and N is a large positive integer 
number). The limit cycle model is obtained by adding an imaginary part to the same 
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asymptotically free model Hamiltonian. The operator T can be replaced by any Hermitian 
operator one wants to use and the methods we employ in this study can still be applied. 



II. DETAILS OF EQUATIONS 



The equations for a fixed value of / contain the diagonal matrix Gf(s) given in Eq. ^ 
that contains the diagonal part of Hf(s). The Hamiltonian matrix is split into its diagonal 
and off-diagonal parts at every value of s, 



H f (s) = D f (s) + V f (s). (10) 



This splitting implies also 



H f (s) = G f (s) + [H f (s)-G f (s)) (11) 
= Gf(s) + f[Df(s) —T] + Vf(s) . (12) 

This means that the diagonal part of the interaction is included in Df(s). The important 
point is that only Vf(s) has non-zero off-diagonal matrix elements. Diagonal matrix ele- 
ments of Vf(s) are zero. T and Df(s) are diagonal and we use an abbreviated notation 
for their matrix elements: T mn = T m S mn , D fmn (s) = D m 5 mn , and G fmn (s) = G m 5 mn . Our 
abbreviated notation for interaction matrix elements is Vf mn (s) = (1 — 5 mn )V mn , where V is 
the full interaction part in the matrix Hf(s) = T + V. 
The SRG Eq. Q implies 

— D n = 2 y]((r n — Gk)V n kVkn ? (13) 

-j^V m7 L n = — (G m — G n )(D m — D n )V mn + (G m + G n — 2Gk)V m kVk n . (14) 

We solve these equations numerically, starting from 

H f (0) mn = y/ E m E n [6 mn - g-ih sgn(m - n) ) , (15) 

E n = b n , (16) 

M < n < N . (17) 
In the case of the model with asymptotic freedom, the coupling constant h is set to 0, 



the Hamiltonian matrix is real and Eqs. (13) and (14) display all relevant formulae. In 



the case of the model with a limit cycle, the Hamiltonian contains an imaginary part and 
the equations we use to compute the Hamiltonian flow numerically involve the imaginary 
components, 

where r is a real symmetric matrix and c is a real antisymmetric matrix. Thus, D m = r mm . 
With this notation, the SRG equations we solve are 

±D n = 2j2(G n -G k )(rl k + cl k ), (19) 

k 

f'm^n {G m G n )( y D m D ri )T mn -\- ^ ^ (G rn + G n 2Gj.) {Trnk^nk ~t~ CnkCmk) i 

(20) 

Cmn (Cm G n j(^D m D n ^jc mn -\- ^ ^ {G rn ~\- G n 2G^,)(c Tra ^r ra ^, c n fcr m /j) . 



ds 



(21) 



We will show that to a good approximation Hf(s) retains the form of Hf(0) for matrix 
elements between states with kinetic energies well below the SRG parameter A = 1/ y^i. The 
dominant change that occurs in these matrix elements is that g is replaced by gf(s). This 
feature emerges in numerical calculations. We also analytically derive this result in the limit 
of large b in the next section. 

Once it is established through numerically calculated non-perturbative evolution of the 
whole Hamiltonian that its evolution can be reduced to the evolution of the coupling constant 
gf{s), we focus discussion on the evolution of the coupling and its dependence on /. The 
coupling constant is defined as 

g f (s) = l-H f (s) MM /E M , (22) 

where Em is the smallest energy in the matrix representation of the theory (smallest allowed 
eigenvalue of T) and Hf(s)MM is the smallest energy diagonal matrix element of Hf(s), the 
infrared corner of the Hamiltonian. We can choose any diagonal or off-diagonal matrix 
element near the infrared corner of the Hamiltonian matrix to define the same coupling 
constant and no significant changes result until this element begins to freeze at its asymptotic 
limit as s — > oo. Each matrix element follows a universal trajectory until it freezes at its 
s — > oo value. Only the lowest diagonal element displays the full evolution of the diagonal. 
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When G differs from D, i.e., / > 0, we show that the SRG transformation does not 
necessarily bring the Hamiltonian matrix to diagonal form as s —>■ oo. However, there is a 
sufficient condition for cases with / > that will be derived in Section [TV] Namely, 

where the derivative means the rate of change of matrix elements of V along the diagonal 
in units of rate of change of eigenvalues of T along the diagonal. This condition could be 
violated when bound-state (negative) eigenvalues appear on the diagonal among positive 
eigenvalues. To study how the SRG transformation behaves depending on the choice of / 
in the presenece of bound states, and in particular what happens in the case / = 1 that 
is useful in nuclear physics, one needs to see what happens in the generic models when / 
deviates from and reaches 1. Since the sufficient condition could be violated when the 
SRG parameter A passes the momentum scale of binding, one should find out how Hf(s) 
behaves around s corresponding to this region. 

In the next section we discuss results of calculations of Hf(s). We show that Wegner's 
transformation (/ = 0) encounters no difficulty and places bound-state eigenvalues on the 
diagonal when A approaches the appropriate bound-state scale from above. The / = 1 trans- 
formation breaks down as the effective cutoff approaches a bound-state scale. In fact, the 
transformation moves the bound-state eigenvalue to the infrared corner of the Hamiltonian. 
When the minimal momentum scale is much less than the bound-state momentum scale, 
the / = 1 transformation forces off-diagonal matrix elements to diverge in order to move 
the bound-state eigenvalue to such a significantly wrong scale. 

For values of / between these two extremes the transformation puts the bound-state 
eigenvalue on the diagonal at some scale between the natural one and the infrared cutoff. 
However, for / approaching 1, the amount of shift approaches the maximal possible and 
correspondingly large couplings must be generated on the way. For / = 1, the transformation 
becomes numerically unstable near the bound-state scale when the infrared cutoff tends to 
zero. 

In discrete matrix notation, one can ask how large a value of / > causes the first shift 
of the bound-state eigenvalue down the diagonal, just by one free energy level in comparison 
to the Wegner case (/ = 0). This value of / will be called critical, and denoted by f c . It 
can be computed numerically to high accuracy. 
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III. NUMERICAL RESULTS: ASYMPTOTIC FREEDOM 



The characteristic SRG behavior displayed by the full matrix Hf(s) is difficult to grasp 
(in our typical cases, the matrix has about 2000 matrix elements, each a function of s that 
ranges from to oo) without studying a number of cases. Such studies involve large amounts 
of data. Fortunately, the overall result of such studies is that the essence of what happens 
for large matrices can be explained in a simple way using just one running coupling constant. 
Most interestingly, the characteristic behavior of this one running coupling constant in large 
matrices can be explained using much smaller matrices, and a matrix that is only 3x3 in 
size will be sufficient. These two facts guide the way we present and discuss our results. 

Our goal is to describe the most important qualitative features of the Hamiltonian evolv- 
ing from its initial form of the type y/E m E n (6 mn — g) to complete diagonalization. We 
need to cover a large range of scales and make important features of V at all scales visible 
simultaneously. Taking our cue from the initial Hamiltonian and experience gathered in ob- 
serving many examples of the SRG flows, we display what we will call the scaled interaction, 
Vmn = Vmn(s) / y/E m E n . Matrix elements of the scaled interaction are all of 0(1) or decay to 
during the entire process of diagonalization when Wegner's transformation is used (/ = 0). 
When / = 1 however, matrix elements of V diverge as the bound-state threshold is reached 
(see below). 

Throughout its evolution, the scaled interaction matrix can be approximated by V mn ~ 
—Qf(s) + corrections for subscripts m and n below the transition region, in which the 
subscripts m and n take values k such that sE% ~ 1. Well above the transition region the 
Hamiltonian matrix is diagonalized. As s increases, gf(s) increases to a maximum, at which 
point the bound-state negative eigenvalue emerges on the diagonal or a process of shifting 
the bound-state eigenvalue down the diagonal begins. 

We start with the model case of asymptotic freedom (a discretized 5-function in two 



dimensions), for which h in eq. ( 15 ) is zero and the single bound-state eigenvalue determines 
g through dimensional transmutation. In the notation of Section [TlJ c mn = and r mn (s = 
0) = VE m E n {5 mn - g). We consider g = 0.0400022797581725654. (This particular value 
of g is not significant; it is found from the condition that for h = 0, b = 4, N = 16, 
M = —25, the bound-state eigenvalue is with high accuracy the same as one of the bound- 
state eigenvalues in another case: one with a limit cycle for g = and h = tan7r/50, and 
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eigenvalue E ~ -7.644479 10~ 6 , see [12].) 




FIG. 1: SRG evolution of H = T + V with A for / = 0. Successive frames correspond to entries 
in Table [ij The horizontal axes display subscripts i and j (running from M to N, counting from 
left to right and from front to back) of the matrix — V and the vertical axis displays corresponding 
matrix elements — V{j divided by ■ s fEjWj. See text for details. 

In Fig. [TJ we display frames from a movie of the evolving interaction, V(s), in the 
Hamiltonian H = T + V, using Wegner's transformation, / = 0. We show —V mn . Since 
small-energy matrix elements of V equal the negative of the coupling constant, we need to 
display —V instead of V itself in order to directly show how well the one coupling constant 
<7/=o(s) approximates the evolution of V below the transition region and at the same time 
show how the coupling constant itself evolves. With these rescaling and display conventions, 
T is the identity matrix. To get H mn j ^E m E n , simply take the negative of the displayed val- 
ues of —V and add peaks of height one along the diagonal in each frame. Details concerning 
the frames shown in Fig. [T] are listed in Table [TJ 

The negative of the scaled interaction is initially a featureless plane (frame 1), its size 
fixed at about 0.04, just above zero because it is a negative of a momentum representation 
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TABLE I: SRG parameters g and A for frames shown in Fig. [TJ numbered from the top to bottom. 
In this example, g = 0.040002, h = 0, b = 4, M = —25, N = 16, and all displayed numbers are 
rounded to 6 decimal places. 



frame 


ln(A)/lnfe 


9(A) 


1 (top) 


22.780321 


0.040002 


2 


2.766096 


0.092055 


3 


-6.864809 


0.600768 


4 


-8.369638 


1.234710 


5 


-8.570281 


0.891475 


6 


-9.071891 


-0.680443 


7 


-12.282193 


-0.226083 


8 (bottom) 


-27.330482 


-0.060769 



of an atractive delta-function in position representation with the initial value of the coupling 
constant g ~ 0.04. As s increases from 0, the plane drops to zero at highest energies, creating 
a cliff (frame 2) between the low and high energy parts of the matrix (transition region). 
This cliff runs along a single row and column that meet on the point along the diagonal at 
which newly decoupled eigenvalues are emerging. The cliff moves towards lower energies and 
grows in height as s increases (frame 3), showing evolution of the transition region between a 
flat plane of zeroed high-energy off-diagonal matrix elements of V and a low-energy plateau 
that is rising higher as s increases (that this plateau rises means that the coupling constant 
increases and the potential itself becomes more negative). The positive peaks left along the 
high-energy diagonal decouple from the rest of the matrix when the ridge moves past and 
they are basically left at the height of the low-energy rising plane at the decoupling point. 

This process continues, with the low-energy part of the matrix rising and high-energy 
eigenvalues being left in isolation on the diagonal as the growing cliff separating low and 
high energy portions of the Hamiltonian moves towards the infrared corner of the matrix 
and the off-diagonal high-energy part of the matrix settles to zero. 

We display a frame in which many high-energy eigenvalues are in place and the low- 
energy plateau has barely fallen from its maximum height (frame 4), at a point where a 
bound state is going to emerge on the diagonal of H(s). Remember that what we are 
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showing is —V, which must be multiplied by \jE m E n and subtracted from T to produce 
H(s). Two important things happen: the low-energy plane is high enough to cancel T and 
produce a negative eigenvalue on the diagonal, and as this bound-state eigenvalue is left on 
the diagonal the low-energy plane reverses its motion (frames 4 and 5) and drops rapidly to 
negative values (frame 6, the coupling constant becomes negative and the whole interaction 
~ gf(s) becomes positive, i.e., repulsive). 

After the bound-state threshold is crossed, V adds to the eigenvalues of T, while up to this 
point it had subtracted from these eigenvalues. After quickly reaching its deepest level, the 
low-energy plane gradually rises to zero (frames 7 and 8), leaving a sequence of eigenvalues 
that smoothly emerge from this flow one after another toward the infrared corner. 

The only violent changes in these scaled variables appear around the point where the 
bound state emerges. This is where the coupling that characterizes the evolution of the low- 
energy part of the matrix grows to its maximum value, before dropping rapidly to negative 
values. 

It is clear that the evolution of the entire matrix V is well-described by the evolution of 
gj(s) in the case of / = 0. The same is true in all other cases we consider. Therefore, we 
will present only the functions gj(s) in all these cases rather than frames from many movies. 

The behavior of coupling constants in the case of asymptotic freedom for various values 
of / is shown in Fig. [2j In the / = 1 case, evolution of H(s) is nearly identical to the / = 
evolution until the bound state begins to emerge. The low-energy plane in V moves up to 
produce the bound-state eigenvalue, but the bound state does not decouple from further 
evolution and the low-energy plane simply continues to grow apparently indefinitely and our 
numerical calculations fail to converge. 

The / = 1 transformation becomes singular because it tries to move the bound-state 
eigenvalue to the lowest-energy diagonal matrix element of the Hamiltonian, as we discuss 
below and explain in greater detail in the next section. This means that the bound-state 
wave function is being forced to include only disparately small momenta. As a result, the 
interactions are forced to grow to maintain observables at their true values. If the evolution 
is halted before the bound-state emerges, the pathological rise of interaction terms does not 
occur yet. 

One can see in Fig. [2] not only that the larger / the larger the maximal value of the 
coupling constant gf at the corresponding value of the argument In A// In b = tf, but also 
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FIG. 2: The coupling constants gj in the case of asymptotic freedom, plotted as a function of 
In A/ In b (instead of s = 1/A 2 ) for 6 values of /: / = (Wegner), / = 0.2, 0.5, 0.75, 0.9, and 
1. The correspondence between a curve and / is such that the curves for larger / reach higher 
and for / = 1 the corresponding curve apparently shoots to infinity around A ~ |-E&<nmdsta*e|- The 
ultraviolet cutoff is at 6 16 , and b = 4. 

that simultaneously the argument tf where the maximum of gf is reached decreases when / 
increases. In fact, when / increases sufficiently to cause a shift of the bound-state eigenvalue 
by one state down the diagonal in comparison to the case / = 0, the maximal coupling 
constant must increase by factor b, when / increases sufficiently to shift the bound-state 
eigenvalue by two states down the diagonal, the maximal coupling constant must increase 
by factor b 2 , and so on. Since Fig. [2] concerns the case with 6 = 4, the maximal values of 
the coupling are about 4 and then 16 when the bound-state eigenvalue is shifted by one and 
two states down the diagonal, respectively. The next shift requires gf ~ 64 and this leads 
to the coupling that grows out of proportion in Fig. [2| When / = 1, the shift occurs to the 
minimal possible value £1 = M, and this requires that g% reaches fr' M ' Eb oun d s tate ~ b 25 ~ 9 = 4 16 
in the case illustrated in Fig. [2] With typical machine limitations, numerical calculations 
are expected to fail below the scale of binding if one uses the generator F{H(s)} with G in 
which / = 1. 
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IV. CONVERGENCE PROBLEMS AND THE 3x3 MATRIX 



In this section, we explain the phenomenon of rise of the coupling constant when the SRG 
parameter approaches the scale of binding in the asymptotically free model. After that we 
will discuss what happens in the model with limit cycle. 

We reduce the asymptotically free model to a 3 x 3 matrix, allowing b to become arbitrarily 
large (i.e., there are only three, strongly coupled degrees of freedom of drastically different 
momentum scales), and show an analytic analysis that explains the behavior observed in the 
full SRG calculation for the 3x3 matrix. The 3x3 matrix model explains the mechanism 
that is also at work in large matrices in our models with asymptotic freedom or limit cycle. 

In the 3x3 case, there is a low-energy effective coupling that evolves smoothly as the 
highest energy scale is decoupled by the SRG transformation. At a characteristic "time" s, 
a 2 x 2 low-energy effective Hamiltonian emerges, characterized by <?/(s), and this coupling 
exceeds some critical value when a bound state emerges. For / = f c , the sign of the dominant 
term driving further off-diagonal evolution changes, and f c can be computed analytically in 
the simple 3x3 example. Above f c the remaining off-diagonal matrix element is forced to 
diverge as a power of b to force the low-energy diagonal to accommodate an eigenvalue of the 
wrong magnitude (wrong in the sense that it is much larger in size than the corresponding 
eigenvalue of T). The b — > oo limit of the full 3x3 matrix evolution can be analyzed 
analytically, and this is how we explain the mechanism at work in the full SRG evolution. 

Before we analyze the drastically simplified 3x3 matrix truncation of the asymptotic 
freedom model, we recall Wegner's demonstration [7j that his transformation always diago- 
nalizes Hamiltonians, and we use similar reasoning to explain why a simpler transformation 
might fail to converge. 

For any similarity transformation, Tr(H 2 ) is independent of s. Separating diagonal and 
off-diagonal contributions to this trace, one finds 



If the magnitudes of diagonal matrix elements increase, the magnitudes of off-diagonal matrix 




(24) 
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elements must decrease. Using eq. (13) for the evolution of the diagonal matrix elements, 

J s E R2 mm = 4 D ^ ~ G n) I Knn \* (25) 

m mn 

= 2j2(D m -D n )(G m -G n )\V mn \ 2 . (26) 

mn 

For Wegner's transformation, G m = D m , so we have 

d 



ds 



Y,H 2 mm = 2 J2(D m - D n f | V mn | 2 . (27) 



No term in this sum can be negative, so the only way Wegner's transformation can stop 
driving off-diagonal matrix elements to zero is if degeneracies appear on the diagonal. In 
this case, the matrix is driven to block diagonal form with diagonal degeneracies in any 
non-diagonalized blocks. 



In general, from Eq. (26 ) we see that negative terms cannot appear on the right- hand- side 



when all differences D m — D n and G m — G n always satisfy the condition 

(G m -G n )(D m -D n )>0. (28) 

Introducing (AT) mn = T m — T n and (AV) mn = V mm — V nn , one obtains for every pair of 
diagonal elements number m and n, 

[/AT + (1 - f)(AT + Ay)] (AT + AV) > , (29) 

which, by dividing by AT > for m > n, implies that v = AV/AT must satisfy the 
condition 

[f + (l-f)(l + v)){l+v)>0. (30) 

This condition implies for / £ [0, 1] that either v < y— j- or v > —1 and only for / = 
(Wegner's generator) these two regions can join while for / > they are always disjoint. 
Instead of differences for arbitrary m and n, it is sufficient to consider differences with 
m = n + 1, since all differences can be built from the differences between neighboring entries 
on the diagonal. Then, in the limit of continuous energy variable, the limit of AT — ► 
produces the condition 

dV 1 dV 

df^T^i or dr-- 1 - (31) 
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as a sufficient one for the transformation to always bring H near the diagonal (outside 
regions of degeneracy mentioned earlier). 

When / = 1, G m = T m , and we see that convergence can fail if 

Ay 



AT 



< -1 



(32) 



for some momenta. T m increases monotonically with m, so we see that problems can appear 
if V m decreases rapidly with m in some region. This is exactly what happens when a negative 
value appears on the diagonal, signaling the appearance of a bound-state threshold. The 
appearance of negative values on the diagonal does not guarantee that the transformation 
will stop driving off-diagonal matrix elements to zero, and it gives no indication that off- 
diagonal matrix elements will actually start to diverge, but it indicates how problems can 
arise. 

To gain further insight, we proceed to a study of the b — > oo limit, which drastically 
simplifies the couplings between various scales, and for further simplicity we truncate the 
asymptotically free Hamiltonian model to a 3 x 3 matrix. We refer to the three remaining 
scales as high, middle and low, and we write the initial Hamiltonian as: 



H(0) 



(b (A 
1 



-9 



( b Vb l\ 

^ 1 7E 



i 

Vb 



(33) 



/ 



where the first matrix is T and the second is V(0). We choose g so that the single negative 
eigenvalue is 0(1) and should appear in the middle of the final diagonalized matrix in the 
case / = 0. These conditions imply that 0.5 < g < 1 and g must be more than 0(\) away 
from either extreme. The three eigenvalues are (1 — g)b + 0(1), (1 — 2g)/(l — g) + 0(l/b) 
and [(l-3g)/(l-2g)]/b + 0(l/b 2 ). 

We define six couplings in the running interaction, 

( -d h b -g h Vb -g„} 
~9hVb -d n 
-9m 



Vf(s) 



\ 



91 

Vb 



9l 

Vb 



(34) 



/ 



where the diagonal couplings dh, d m and di, and the off-diagonal couplings gh, g m and g\ are 
all functions of s. At s = all of these couplings are equal to g. The scaled interaction V is 
obtained from Vf(s) by replacement of b by 1. 



15 



Our goal is to obtain accurate estimates for all matrix elements of Vf(s), even though 
their magnitudes span many orders. We cannot allow small errors in large eigenvalues or 
far off-diagonal matrix elements to produce large errors in small diagonal matrix elements 
that should reproduce the eigenvalues of order 1 or 1/6, the prototypical renormalization 
problem. As b — > oo this problem can be analyzed analytically. 



We find evolution of the coupling constants in Eq. (34) using Eqs. (13) and (14). The 



full set of equations for arbitrary / is (dots indicate derivatives with respect to s) 

d h = -2b'y hm gfa-2'y hi gft n , 

d m = 2b 2 i hm g\ - 2b~ 1 i hl gf , 

d\ = 2b 2 i M g m + 2<y ml gf , 

9h = ~b 2 7hm[(l - d h ) - b~ x {l - d m )] g h - (ihl + k'Smi) 9m9l , 

9m = ~b 2 7w[(l - d h ) - b~ 2 (l - di)]g m - b(j hm - b" 1 ^) 9h9i , 

9l = b 2 (jhm + Ihl) 9h9m ~ 7mi[(l ~ d m) ~ & _1 (1 - d i)] g { , 



(35) 
(36) 
(37) 
(38) 
(39) 
(40) 



where 



7/itj 1h 1m/b i 

1ml = 7m - ll/b, 

Ihi = Ih ~ li/b 2 , 

lh = / + (i-/)(i-4), 

7m = f + (i-f)(i-d m ) 

n = / + (l - - d,) . 



(41) 
(42) 
(43) 
(44) 
(45) 
(46) 



A. Approximate evolution of 3 x 3 matrix for / = 

We begin by considering Wegner's transformation, / = 0. We use Eqs. (|35|)-([40|) and 
keep only the leading terms for large b. We will see that all couplings remain 0(1) when 
/ = 0, so this analysis is fairly straightforward. The evolution has two stages: elimination 
of gh and g m in the first stage, and elimination of gi in the second stage. Terms driving 
the first stage of evolution are 0(b 2 ) and govern evolution until s exceeds 0(l/b 2 ). These 
terms are then exponentially suppressed by gh and g m , and the second stage of evolution is 
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governed by subleading terms. The leading terms are: 

d 



ds 
d 



d h 



~dr, 



ds 
d 

ds 
d 

ds 
d 



di 



9h 



ds 



9u 



d 

ds 



9i 



— A{i — a h )g h o + , 


I 4 'J 


rl, \n 2 h 2 -1- D(h\ 




2(1 - d h )g 2 m b 2 + 0(1) , 


(49) 


-(l-d h ) 2 g h b 2 + 0(b) , 


(50) 


-(l-d h ) 2 g m b 2 + 0(b) , 


(51) 


2(1 - d h )g h g m b 2 + 0{b) . 


(52) 



We see from (50) and (51) that g^ = g m during the first stage, because all couplings start 
at g. This implies that d m = di = gi during this stage also, because the leading equations 
governing their evolution become identical. Since the low energy 2x2 submatrix of V is 
determined by d m , d\ and gi, it retains its original form during the first stage of evolution, 
with a single coupling that can be factored from the submatrix of V. This is one of the most 
important results of this analysis and it can be generalized to larger matrices. 

The leading term on the right of Eq. (47) is 0(b) rather than 0(b 2 ), so for s of 0(l/b 2 ), 



dh changes only by 0(l/b) and we can ignore this change when solving Eqs. (50) and (51), 



replacing dh with its initial value, g. Solving Eqs. (|50|) and (|51|) for large b we obtain 

9h(s) 



9m{s) ~ g e ^ a ' - " . (53) 
Both of these couplings decay to zero exponentially, and for large b this decay is so rapid 



that it can be treated as instantaneous. Inserting Eq. (53) in Eq. (47) we find 

dk = 9~ 7V~~XZ + 0{l/b 2 ) , 



(54) 



at the end of the first stage of evolution. For a complete leading-order analysis, we need 
only these leading approximations for gh, g m and dh- The largest eigenvalue, (1 — dh)b ~ 
(1 — g)b + 0(1) thus appears in the high-energy corner of the matrix. We do not need the 
0(l/b) correction to dh to obtain the highest eigenvalue accurately for large b and we will 
see that we do not need this correction to accurately compute the smaller eigenvalues either. 
The leading term governing the early evolution of the remaining 2x2 submatrix is 

2 



2(1 



d h )g 2 h b 2 ^ 2(1 - g)g 2 b 2 e- 2 ^ 2b2 ° 
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9 



9 



-As) 



(55) 



where S(s) is defined so that J °° S(s) ds — 1. The couplings d m , di and gi instantly increase 
by g 2 /(l — g), so at the end of the first stage and the beginning of the second stage of 
evolution 

d m {s) » d t (s) » » ^- + 0(1/6) . (56) 

For the range of (/ that produces a binding energy of 0(1), these low- and middle-energy 
couplings exceed 1 and negative values appear on the diagonal of H at the start of the 
second stage of evolution. 

After the initial instant of evolution, g h « 0, g m « and dh ~ <?■ Corrections to these 
approximations have no effect on the leading order of any eigenvalue. For example, the 
0(l/b) correction to dh has no effect on the smallest eigenvalue, which is 0(l/b). Returning 
to the full equations, replacing g h , g m and dh with these approximations, choosing an initial 
value of g/(l — g) for d m , di and gi, we find the leading-order equations that govern the 
second stage of evolution: 

^-d m w -2(1 - d m )gf l - « , (57) 
^d, « 2(1 - d m )gf , (58) 
i. gi » _(i_ dm )2^ . (59) 

When (7, the initial coupling, is chosen so that the binding energy is 0(1), 1 — d m < at 
the start of the second stage and for large b, d m does not evolve further. This means that 
the negative eigenvalue, 1 — d m — 1 — g/(l — g) — (1 — 2g)/ (1 — g) appears in the middle of 
the matrix and stays there. This leaves us with, 

Ui - 2 ^gf , (60) 
ds 1 — g 

d ( 1 2g s : 

From these we see that 



ds 91 * " lT^7 1 ^ • (61) 



9l (s) » e-^-^VCi-s) 2 s , (62) 



which leads as s — > cxo to 



This means that the smallest eigenvalue, (1 — d{)/b w [(1 — 3g)/(l — 2g)]/b, appears in the 
infrared corner of the matrix when / = and no couplings become unnaturally large during 
the evolution. 
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B. Approximate evolution of 3 x 3 matrix for / = 1 



Next we study the transformation when / = 1. We again use Eqs. (35)-(40) and keep 
only the leading terms for large 6 for the first stage of the evolution. We will see that and 
g m again decay exponentially and all other couplings remain 0(1) during this first stage, 
but when / = 1 the coupling gi grows during the second-stage to 0(6 1//2 ) before it is finally 
driven to zero and the second-stage's analysis is much more complicated because of this 
unnatural growth. The leading terms for / = 1 are: 

±d h = -2^6 + 0(1), (64) 

j- s d m = 2 9 y + 0(6) , (65) 

j- s d t = 2g 2 J? + 0(1) , (66) 

^g h = -(l-d h )g h b 2 + 0(b), (67) 

^g m = -(l-d h )g m b 2 + 0{b), (68) 

^-9i = 2g h9m b 2 + 0(6) . (69) 
as 

These equations govern the first stage of the SRG evolution, with subleading terms becoming 
important in the second stage. We again see that g^ = g m to leading order, which implies 
that d m = di = gi to leading order during the first stage. The low energy 2x2 submatrix 
is determined by these three couplings, and it retains the form of the initial Hamiltonian 
for all values of / during the first stage of evolution, so we can use the universal coupling 
governing the submatrix to characterize the transformation for all values of /. Solving the 
equations for g h and g m for large b, using the fact that dh = 9 + 0(1/6), we obtain 

9h(s) « g m (s) « 9 e-M 62 * . (70) 



This is nearly the same result we obtained in Eq. (53) for / = but the exponent is 
different. Both couplings decay to zero exponentially and dh changes only at 0(1/6) because 
the exponent is 0(6 2 ) while the term driving the evolution of dh is only 0(6). The largest 
eigenvalue, (1 — dh)b w (1 — g)b again appears in the high-energy corner of the matrix. 
The leading term governing the early evolution of the remaining 2x2 submatrix is 

2^6 2 » 2« 7 2 6V 2 ( 1 -^ » ^-S{s) , (71) 
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where 8(s) is again defined so that / °° 8(s) ds = 1. Once again we find that the couplings 
d m , d[ and gi instantly increase by g 2 /(l — g), to leading order, so at the beginning of the 
second stage of evolution 

d m (s) « « <?z(s) « ^- . (72) 

These second-stage "initial" values can differ at (9(1/6), but such differences do not affect 
the subsequent analysis when / = 0. We will see that a full analysis when / = 1 is sensitive 
to such 0(l/b) corrections, but we are not interested in the effects of these small corrections. 
We focus on how the transformation works when the simple scaling analysis breaks down. 
We want to understand the origin of this breakdown and we can infer its consequences 
without a derivation of the precise evolution of the couplings. The origin lies in the specific 
SRG transformation, not in the precise second-stage initial values of the coupling constants. 

For the range of g that produces an 0(1) binding energy, the low- and middle-energy 
couplings exceed one and negative values again appear on the diagonal of H. However, when 
/ = 1 the appearance of negative values on the diagonal can signal trouble, as argued at the 
beginning of this section. 

After the initial stage of evolution, « 0, g m ~ and dh ~ g. The subleading terms 
that govern the second stage of evolution in this case are: 

jd m « , (73) 

£«■ « 2 9? , (74) 

d , 1 — d\ . . 

—gi « -(1 - d m )gi H —g l . (75) 



When the binding energy is 0(1), at the start of the second stage 1 — d m < in Eq. (75), so 



gi grows exponentially rather than decaying as it did when / = 0. Eq. (73) implies that d. 



will decrease monotonically at 0(1/6), but before it can decrease sufficiently to reverse the 



growth of gi, g\ grows to 0(\/b). Meanwhile, Eq. (74) implies that di will grow monotonically 
and it eventually grows to 0(6) as the bound state eigenvalue is moved from the middle of 
H to its low-energy corner. Once d m decreases sufficiently and d\ grows sufficiently, the sign 



of the right-hand-side of Eq. ( 75 ) changes and gi is driven to zero exponentially from this 
point. 

An exact solution of these equations is readily found because the unitary evolution of the 
2x2 low-energy submatrix of H is simply a rotation that conserves its trace and determinant. 
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Since d m decreases monotonically, the bound state eigenvalue cannot appear in the middle 
of the matrix, as it did when / = 0. The low-energy eigenvalue that is (9(1/6) must end up 
in the middle of the matrix, and this means the precise final result of the transformation 
is sensitive to 0(l/b) corrections to the initial value of d m , but the mechanism by which 
the 0(1) eigenvalue is transferred to the 0(l/b) momentum corner of H is not. Thus, the 
unitarity of the transformation combined with the orders of magnitude of the couplings 
d m , gi and di at the beginning of the second stage are enough to diagnose why numerical 
transformations of large matrices fail to converge for some transformations. Moving the 
bound state eigenvalue to an "unnatural" location requires growth of couplings by powers of 
b, and in large matrices additional powers of b appear each time the bound state eigenvalue 
is moved one step down along the diagonal. 



C. The critical value f c 



We do not provide complete details of the full evolution for arbitrary / between zero 
and one; the critical value, f c , is revealed at the end of the first stage of evolution. As the 
bound-state eigenvalue emerges on the diagonal, / determines whether gi grows or decays 
exponentially. 



We again use Eqs. (35)-(40) and keep only the leading terms for large b. 



d 

ds 
d 



d h 



(I,,. 



ds 
d 

ds 
d 

ds 
d 



-di 



9h 



ds 



9n 



ds 



9i 



-2[1 


-(l-f)d h ]glb + 0(l), 


(76) 


2[1- 


(l-f)d h ]g 2 h b 2 + 0(b), 


(77) 


2[1- 


(l-f)d h ]g 2 m b 2 + 0(l), 


(78) 


"(1- 


-d h ) [l-(l-f)d h ]g h b 2 + 0(b) , 


(79) 


-(1- 


-d h ) [l-(l-f)d h ]g m b 2 + 0(b) , 


(80) 


2[1- 


(l-f)d h ]g h g m b 2 + 0(b). 


(81) 



As expected from the two extremes / = and / = 1, during the first stage of evolution, 
9h = 9m, implying that d m = di = g%. Moreover, the entire analysis of the first stage of 
evolution is qualitatively independent of /. To leading order: 



9h{s) = g rn {s) = g e-MMW)^ 



(82) 
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cf. Eqs. g) and 

Once again the integrated effect of these couplings on the low-energy 2x2 submatrix is 
to shift the submatrix couplings to a new starting value, at s = + = 0(l/6 2 ), 

(7,(0+) = d m (0+) = d,(0+) = ^- , (83) 

after which we need only examine the equation governing subsequent evolution of gf. 

^gi = - [1 - d m (s)} [1 - (1 - f)d m {s)} 9l {s) + 0(1/6) . (84) 

This off-diagonal coupling grows for d m {s) > 1 if 

l-(l-/) I ^->0, (85) 

from which we determine 

fc = 2 - - ■ (86) 

This value is always between zero and one if there is a bound state with a binding energy 
that is 0(1). 

The evolution of larger matrices as b — ► 1 is far more complicated than the b — ► oo limit 
of the 3x3 matrix, but there is remarkable similarity. Any stage in the evolution of larger 
matrices can be modeled by a 3 x 3 matrix with the middle diagonal matched to the point 
on the diagonal where an eigenvalue is emerging. We have focused on the point where a 
bound-state eigenvalue appears, because in the evolution of large matrices with b near 1, it is 
at this point in the SRG evolution of the Hamiltonian that transformations bifurcate; those 
with f < f c leave the bound state eigenvalue on the diagonal where it appears for / = 0, 
while those with f > f c move it down the diagonal. Transformations with / — > 1 misplace 
the bound-state eigenvalue to such small scales that they fail to numerically converge. Off- 
diagonal matrix elements are forced to diverge exponentially to accomplish this. There 
is always a tipping point, f c < 1, and at f c when the bound state emerges, off-diagonal 
matrix elements are balanced at an 0(1) value between regions of exponential decay and 
exponential growth. 



V. NUMERICAL RESULTS: LIMIT CYCLE 



The limit cycle Hamiltonian of Eq. (15) with h ^ possesses many bound states, one for 
each cycle. The number of eigenvalues that emerge in each cycle is fixed by h. Application 
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of the SRG procedure in the case of a limit cycle is based on [12]. Numerical calculation 
in this case produces Fig. [3} One sees that the case of asymptotic freedom in the previous 
subsection corresponds in its behavior around the scale of binding to one cycle around 
the scale of binding in the limit cycle. In addition, one sees that when f —> 1 and the 
SRG generator becomes purely kinematic, the RG evolution is stuck in the range of scales 
corresponding to the scale of binding in the first cycle (greatest binding energy). In order 
to get to the next cycle and smaller binding energies, one has to introduce dynamics into 
the generator [G, H] through the operator G. 

Note that Fig. [3] also shows that the larger the shifts of the bound-state eigenvalues 
down the diagonal and the larger the corresponding coupling constant, the more sensitive 
the numerics are to details of the finite matrix, introducing departures from a clean cycle. 
These numerical effects are not fully understood but they have no bearing on the findings 
reported here that concern the impact of bound states on the SRG transformations. 

It is now clear that in the cases where the greatest binding energy E is very small, one 
does not encounter problems using purely kinematic G = T as long as A ^> \E\. This is our 
explanation of what happens in the studies in nuclear and atomic physics where the binding 
energies are much smaller than the cutoffs on momenta. Hadronic structure is different 
because there are no free quarks or gluons; all high energy eigenstates involve bound states 
of quarks and gluons so it is not possible to remain above a bound-state threshold. 

When one proceeds to the case of a limit cycle with a geometric series of binding energies, 
one must insure that the largest of them is small in comparison with desired values of A for 
the generator with G = T to work. 

VI. CONCLUSION 

The SRG offers an alternative to traditional renormalization group transformations which 
discard degrees of freedom and is being developed to attack problems where traditional 
methods fail. These problems all involve a broad (in principle, considered infinite) range of 
momentum scales that are strongly coupled and are typically not amenable to perturbation 
theory, even if the latter is so-called renormalization-group improved. Analytic methods 
fail and direct numerical solutions are not possible because the space of states cannot be 
truncated, due to the strong coupling, and the full space is too large for current or foreseeable 
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numerical storage and manipulation. In the most interesting problems one cannot even 
be sure what operators are needed to construct renormalizable Hamiltonians (e.g., what 
counterterms are required in the presence of known bare interactions). 

We have not illustrated the initial steps required to find renormalizable Hamiltonians, 
choosing simple model problems in which the required operators are basically known (e.g., 
see [23]). But our calculations do illustrate methods that can make such initial steps feasible. 
First, as anticipated in Wilson's earliest work [24], we drastically truncate the space of states 
by moving from a continuum to a space of discrete states which are spaced exponentially, 
with a tunable spacing governed by a parameter b. Going from a well-defined continuum 
problem to a discretized problem is straightforward in the models we use, so we have paid 
little attention to this step; nor have we shown that with the renormalization problems we 
solve the continuum can be recovered with exponential convergence by letting b — > 1. The 
relevance of our calculations rests on the assumption that such a limit can be taken with 
control of errors and without introducing new renormalization problems. Once we identify 




FIG. 3: The coupling constants 5/(A) in the case of a limit cycle, plotted as a function of In A/ In b 
for 6 values of /: / = (Wegner), / = 0.2, 0.5, 0.75, 0.9, and 1. The correspondence between 
a curve and / is such that the curves for larger / reach higher and for / = 1 the corresponding 
curve apparently shoots to infinity already around In A/ In b ~ 15. The ultraviolet cutoff is at 6 16 
and 6 = 4, see [T2] . 
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the operators required for renormalization with large 5, using nonperturbative numerical 
calculations like those shown here, we reach what is the starting point for the calculations 
we present. 

Attempts to control the effects of strong coupling over large numbers of momentum scales 
typically fail because divergences emerge. Originally, such divergences showed up in pertur- 
bation theory and gave rise to the whole renormalization program of the last century. Far 
more challenging are the types of divergence illustrated by our calculations, divergences that 
persist if one attempts to go beyond perturbation theory using perturbative renormalization 
schemes such as standard methods for the renormalization of quantum electrodynamics. 
Wilson's renormalization group improved perturbation theory avoids many of these prob- 
lems in asymptotically free theories as long as cutoffs are kept sufficiently large. However, to 
solve problems in which strong-coupling fixed points [3] or limit cycles exist, renormalization 
group transformations often must be crafted on a case by case basis. What the SRG offers is 
a wealth of new transformations and in this paper we have focussed on some of the critical 
features of these new transformations. 

The existence of a critical value / = f c in SRG transformations with generators F{H} = 
[G, H] and G = fT + (1 — f)D, is demonstrated in discrete models with important features 
such as asymptotic freedom or a limit cycle. The case of SRG with / = 1, or G — T, is 
always in the region / > f c . Therefore, when bound states exist and the SRG parameter 
A = 1/ \/s approaches the scale of momenta that dominate in the formation of bound states, 
the strength of renormalized interactions grows. Numerical calculations of the interaction 
Hamiltonians become increasingly difficult due to this growth. But if it is enough to calculate 
effective theories with A much larger than the scale of binding, the generator with / = 1 can 
be employed without encountering an intractable increase of interaction strength. 

As long as one stays away from the bound-state momentum scale, the numerical calcu- 
lations in the generic models are equally powerful for / > f c as they are for / < f c . This 
result shows that limits on applicability of the simplest version of the SRG transformation 
with / = 1 due to bound states are not as severe as one might expect provided one keeps 
A away from the scale of binding. This is important because the SRG transformations with 
/ = 1 are the simplest to implement in the continuum limit and in perturbative evaluation 
of the SRG flow of Hamiltonians. 

At the same time, it is also made clear that in order to handle cases with a limit cycle, 
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one has to consider / < f c . This means that the generator must include interactions in 
G and cannot be limited to T. The / = transformation advocated by Wegner [7J [S] is 
able to drive both of our models to diagonal form. We have not shown that this reproduces 
both the correct binding energy and phase shifts for the continuum problem but such a 
demonstration should be straightforward and is left for future work. 

Finally, it should be pointed out that this article does not resolve many important issues 
that must be resolved to deal with confinement. These go well beyond dealing with bound 
states and are not encountered in our simple generic models, but confinement presents us 
with the problem of bound states at all cutoffs and at least in this respect the limit cycle 
model should provide important insights. 
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